Penguin pollution in the B° — > J/ij)Ks decay 
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We present the most complete analysis of the penguin correction to the extraction of the standard- 
model parameter sin(2<^>i) from the P>° — > J/ijiKs decay up to leading power in 1/mj, and to next- 
to- leading order in a s , 4>i being the weak phase, nib the b quark mass, and a s the strong coupling 
constant. The deviation ASj/^ Ks of the mixing- induced CP asymmetry from sin(2^i) and the 

■ direct CP asymmetry Aj^ Ks axe both found to be of 0(1CP 3 ) in a formalism that combines 
the QCD-improved factorization and perturbative QCD approaches. The above results, different 

<3 from those of O(10" 4 ) and of O(10" 2 ) obtained in the previous calculations, provide an important 

■ standard-model reference for verifying new physics from the B° — » J/ipKs data. 

H ■ 

PACS numbers: 13.25.Hw, 12.38.Bx, ll.10.Hi 

The B° — > J/ifjKs decay has been regarded as the golden mode for extracting the standard-model parameter 
sin(20i) [H, cf)i being the weak phase of the Cabbibo-Kobayashi-Maskawa (CKM) matrix element Vtd- Though 
£C) ■ it is believed that the penguin pollution in this mode is negligible, a complete and reliable estimate of its effect 
| is not yet available. Such an estimate is essential, especially when B physics is entering the era of precision 
measurement. Note that the data of the mixing-induced CP asymmetries in B meson decays, such as of Sj/^^s > 
are not included in the CKM or UT Q fitter. Understanding the penguin effect, one will have a more concrete 
idea of whether a discrepancy between the measured Sj/^Ks an d sin(2</>i) from the fitters is a signal of new 
physics. At the same time, the direct CP asymmetries of the B — > J/ipK decays, also related to the penguin 
correction, are expected to be vanishingly small. A reliable estimate of the penguin effect, when confronted 
with the precision measurement of Aqp{B — > J/tpK), can reveal a new physics signal. It has been claimed that 
Acp(B — > J/ipK) observed at 1% or higher would indicate new physics definitely Q]. 

The previous calculation of the deviation ASj/-^k s = Sj/ipK s — sin(20i) ~ O(10~ 4 ) [H has taken into account 
the corrections to the B-B mixing [fjj and to the decay amplitude. However, only the tt-quark-loop contribution 
from the tree operators was included in the latter. We shall point out that the contribution from the penguin 
operators has been overlooked in [f|, which may result in a larger deviation. The model-independent analysis 
in Q gave ASj/^k s = 0-000 ± 0.012, whose theoretical uncertainty is comparable to the current systematical 
error in the Sj/^Ks measurement. This method relies on the input of the B° — > J/tpir data in order to fix the 
penguin amplitude (including its strong phase) . The experimental error then propagates into their uncertainty 
of O(10~ 2 ), which should have been overestimated. Though the B° — > J/ipir data will get precise in the future, 
SU(3) symmetry breaking effects are also going to emerge, and one will lose the control of the theoretical 
uncertainty eventually, if using this method. Hence, we are motivated to reinvestigate this important subject 
by computing the complete penguin correction to the B° — > J/tpKs amplitude. 

In this work we shall adopt a special formalism [8] that combines the QCD-improved factorization (QCDF) 
[§| and perturbative QCD (PQCD) [l(| H3 approaches. The factorizable contribution in the B — > J/tpK 
decays contains the B — > K form factor F+ (m^) evaluated at the invariant mass squared of the J/tp 

meson. Because of the heavy mj/^, the energy release involved in F+ K (rn 2 j^) is small, and it is unlikely to 
further factorize it into a convolution of a hard kernel with the B meson and kaon distribution amplitudes as 
in PQCD. Therefore, we adopt the QCDF formalism to handle the factorizable amplitude. The nonfactorizable 
spectator contribution is essential in the color-suppressed category of B meson decays [TH, [l3|, the reason 
the naive factorization assumption [l4| does not apply to the B — > J/tpK decays. For this piece, QCDF is not 
appropriate due to the end-point singularity from vanishing parton momenta, when the twits-3 kaon distribution 
amplitudes are included [15j . Therefore, we employ the PQCD approach based on kr factorization theorem, 



X 



*Electronic address: hnli@phys. sinica. edu.tw 
^Electronic address: mishima@ias.edu 



2 



which is free of the end-point singularity. As argued in , the nonfactorizable contribution has a characteristic 
hard scale higher than that in the B meson transition form factor. This mixed formalism has been applied to 
B meson decays into various charmonium states [8| , and results for the branching ratios are consistent with the 
observed values. 

The mixing- induced CP asymmetry Sj/^k s anc ^ the direct CP asymmetry Aj^ Ks of the B° — > J/ipKs 
decay are defined by 

2ImXj/^ Ks \^j/j,K s \ 2 ~ 1 

J/ ^ KS ~ 1 + \\ JH ,Ks I 2 ' 1 + \^J/ipK s I 2 ' 

_ q A{B Q ^J/^K s ) 
J/ ^ Ks pA(B° -> J/I/Ks)' [ > 

where the ratio q/p is related to the non-diagonal elements of the mixing matrix, and A the decay amplitude. 
Note that a tiny Kl admixture could exist in the experimentally reconstructed Kg final state, such that CP 
violation in the K-K mixing contributes to Sj/^Ks anc ^ ^J/i>K s - However, we shall concentrate only on the 
corrections from the B meson system in this work, and refer the inclusion of the corrections from the kaon 
system to [l7| . 

We aim at an accuracy up to leading power in 1 /m b and to next-to- leading order in a s , m b being the b quark 
mass and a s the strong coupling constant. At this level of accuracy, the correction to the mixing factor q/p, 
computed reliably in the framework of operator product expansion with the large logarithm a s ln(mw/mb) 
being resummed [5(, gives 

= (2.08 ± 1.23) xlO- 4 , (2) 

where my/ is the W boson mass. Note that the old value sin(20i) » 0.736 has been input to derive the above 
result. Using the updated one sin(20i) w 0.684 [18[ changes Eq. by at most 10%, namely, by O(10 -5 ), which 
is absolutely negligible. The influence of the B-B mixing on the tagging leads to the direct CP asymmetry 

= (2.59 ± 1.48) xHT 4 , (3) 

which is consistent with the estimate in [2(|. Both Eqs. © and (J3]) are of O(10~ 4 ). 
The B Q — ► J/ipK° decay amplitude is decomposed into 

A(B° - J/yjK ) = V;„V CS (A { j% K0 + Af^ K0 ) + v: b V us (A% K0 + A? /ilK0 ) , (4) 

where the unitarity relation V t * b Vt s — —V* b V us — V* b V cs for the CKM matrix elements has been inserted, and the 
second term carries the weak phase <p 3 . The dominant tree amplitude A^ J p K0 is responsible for the B° — * J/i/jK a 

branching ratio, and A!j$ k0 are the penguin pollution in the decay. The amplitude Ajj^ K0 , receiving the 
contribution from the u-quark loop through the BSS mechanism [2lj . has been estimated naively in [f|, while 
Aj,^ K o from the penguin operators has been missed. We shall calculate these amplitudes in our formalism 
below. 

The leading-power QCDF formulas for the factorizable contributions to A^j'^ k0 are written as Q 



Af/^ K0 = 2v/2iVc/ Q dx 3 ^ L (x 3 )a 2 (x 3 ,t)Ff K (m 2 j / ^ , 



Af/^ K0 = 2^/2N c ( dx 3 * L (x 3 )[a 3 (x 3 ,t) + a 5 (x 3 ,t)] F^ K {m 2 j^) , (5) 



with iV c = 3 being the number of colors, ^f L the twist-2 J/tp meson distribution amplitude, and x 3 the momen- 
tum fraction carried by the c quark. The effective Wilson coefficients including the 0(a s ) vertex corrections 
are given by 



a 2 {x,p) = Ci(/i) + ^M 



1 + (-18 + 12 In H» + fl (x) 

47T V U 
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a 3 (x :f i) = C 3 (/i)+C 9 (/i) + — [Cii/i) + Cho(jm)] 
a 6 (x,ii) = C 6 ^)+C 7 ( l i) + — [C 6 (ji) + C a (jj)] 



1 



a s (/x) 



47T 



-18+ 12 In— + fi(x] 



6 -12 In — -/j(x) 
A 4 



where the explicit expression for the function // can be found in [la. |22| | , 

2t|(1-x) 



fj{x) = - ^ In a - 2,tti + 31n(l - r|) 



1 — TgX 



(G) 



(7) 



with the ratio r 3 — vnji^/mB, rnB being the B meson mass. Compared to jg], the hard scale t for the strong 
coupling constant and for the Wilson coefficients has been chosen as the larger internal momentum transfers 
between (P 2 + x 3 P 3 ) 2 w x 3 (m 2 B 



l J/i>- 



and (P 2 + X3P3) w X3(m^ — rriju) with the notation X3 = 1 — £3, 



i.e., 



i = max 



^x 3 {m 2 B - rnj^), ^Jx 3 (m 2 B - m 2 ^) j . 



(8) 



This choice is more consistent with the derivation of the PQCD formalism [23[ and follows the well-known 
Brodsky-Lepage-Mackenzie procedure [24| . 

The leading-power 0(a s ) nonfactorizable spectator amplitudes in PQCD are quoted from 



At f K<> = Mfa) , A^ K0 = M(a> 3 - a' 5 ) , 
where the effective Wilson coefficients are defined by 

40«) = ^ , = ^- [G*0*) + CM] , a' 5 ( M ) = i [q>( M ) + c 8 ( M )] 

The amplitude .M is written as 



,(*)«/ _ 



(9) 



(10) 



16irC F ^/2N c 



dxidx 2 dx 3 







6id6i(/> B (xi,6i) 
1 



x{ [(1 - 2r 2 3 )x 3 -<f ) £(x5)* L (x 3 ) + ^r 2 ^(a^)**(x 3 ) 



-r 2 (l - r 2 )x 2 ^(aJ2-)* L (x 3 ) + r 2 (2rf&3 + (1 - r 2 )^) ^(a^)*^) 



- \x 3 + (1 - 2r 3 2 )x 2 ) ^(x 2 -)* i (x 3 ) + r\ (2r 2 <^(xi) - -^(35))** (a*) 
-r 2 (l - r 2 )x 2 ^(x^l^) - r 2 (2r 2 x 3 + (1 - r 2 )^) ^(si)* 2 ^) 
xa s {tf)a[{tf)S{tf)hf{x u x 2 ,x 3 M)}, 



(11) 



with the color factor Cf = 4/3, the mass ratio r 2 = uiqk /tub, nrioK being the chiral enhancement scale 
associated with the kaon, the impact parameter b\ conjugate to the transverse momentum of the spectator 
quark, the two-parton twist-3 J/ ifj meson distribution amplitude , and the notation xj = 1 — x 2 . The explicit 

expressions of the hard functions h} , the hard scales t^' 2 \ and the Sudakov factor S are referred to 0]. 

The penguin correction from the M-quark loop, which is of next-to-leading order in a s , can be expressed in 
terms of the effective Hamiltonian 



nj(u) _ GF t/ * t/ 

2 , I 2 
In — k 



^e u e c (N e d + C 2 ) {cc)v{sb) V -A + ^C 2 (cT a c) v (sT a b) v _ A 



/7T 



(12) 
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The first (second) operator arises from the process, where the c-quark pair is produced through the photon 
(gluon) emission from the M-quark loop. I 2 is the invariant mass squared of the photon or gluon. Notice the 
constant 2/3 different from 5/3 in Q, since we have considered the Fierz transformation of the four-fermion 
operators with the anti-commuting 75 in D dimensions under the NDR scheme [25j . We have also included the 
appropriate color factor N c and the charge factor e u e c = 4/9. It is easy to observe that the first operator can 
be decomposed into O7 and Og, and the second operator into 03-6- 

The leading contribution from the first operator is a factorizable amplitude, 

A { jH S k0 = 2 VlN~cJ o dx^ L {xz)a ew {xz,t)F* K {m 2 J/%l> ) , (13) 
with the effective Wilson coefficient, 

/ r, ^2 \ 

J7T \ 3 /J, 



x,^) = -— e u e c [N c C 1 ( t j,)+ C 2 (ii)][- -In— ^+ in) , (14) 



'•j/ip- We nave examined the average hard scale for the factorizable contribution, 
which is slightly above 3 GeV (the average hard scale for the nonfactorizable contribution is between 1 and 1.5 
GeV). Simply comparing Eq. (fT4| with Eq. ((6| for n 3 GeV, we find that the magnitude of Eq. (fT3|) is less 
than 5% of \Aj,^ k0 \. Including the higher-order terms, such as the nonfactorizable amplitude, the effect from 

the first operator is expected to be around 5%, and can be safely dropped. 

For the second operator in Eq. (|12p to contribute, an additional gluon must attach the c-quark pair in the 
color-octet state. If it is a soft gluon, the corresponding nonperturbative input is the three-parton J/ip meson 



distribution amplitude, defined via the nonlocal matrix element 26] 



(JMP, (* w )\c(-z)g s G^{vz) la c{z)\0) = ifj^mj^P a [p v e*^ - e^P,,] V(v, P ■ z) 

e *(A) . z 

+ifj/1,m*j /i> [P^gt - P,5a M ] P ■ z) 



P a [P^z v -P v z^{v,P-z) , (15) 



with the polarization vectors e*- A ' of the J/ijj meson, the gluon field strength tensor Ga V , and the projector 
9uu = 9/j.v — (P^Ziy + P lJ z f2 )/(P ■ z), z being the coordinate of the c quark field. Because mj^/ms is not a small 
ratio, we extend our investigation to twist 4. The twist-3 distribution amplitude V , involved only in decays 
into transversely polarized J/tp mesons, is irrelevant here. However, it could be a source of uncertainty for 
the extraction of sin(20i) from the B — > J/tpK* decays. The twist-4 distribution amplitudes $(a; c , x 5 , x g ) and 
f (j; c , x 5 , x g ), from the Fourier transformation of $ and ^, respectivel y, a re anti-symmetric under the exchange 
of the momentum fraction x c of the c quark and x e of the c quark [26j |. x g being the soft gluon momentum 
fraction. Because the hard kernel associated with the u-quark loop is symmetric under the exchange of x c 
and Xc, its convolution with $ and ^ diminishes. That is, the contribution from the second operator through 
higher-power terms is suppressed. 

If the additional gluon is hard, the resultant contribution, being of 0(a 2 s ) (the gluon emitted by the u-quark 
loop is also hard), is beyond the scope of this paper. Nevertheless, we shall have a closer look at this piece. 
When the additional gluon is emitted by the spectator quark, the corresponding amplitude has the formula the 
same as -4j*/^o but with the replacement, 

°>^T£t„ c *(l-'»7 + '') ' (16) 

Since the argument of Aj,^f K0 depends on the difference of 03 and a' 5 as shown in Eq. © , the above corrections 
cancel each other exactly. When the additional gluon is emitted by the b quark, the s quark, or the u-quark 
loop, the corresponding amplitude involves a two-loop calculation, for which there is no simple estimation. The 
two-loop vertex corrections from various operators are expected to be finite, so their complete analysis for both 



5 



^j/ipK ano - ^jjipK ma y ^ e required in the future. Following the above analysis, we shall also ignore the 
c-quark-loop correction, which modifies only the prediction for the branching ratio slightly. 

The choices of the b quark mass, the B meson wave function 4>b, the B meson lifetime, the twist-2 kaon 
distribution amplitude (f> K , the two-parton twist-3 kaon distribution amplitudes <f> K and (f> K , the B meson and 
kaon decay constants, the chiral enhancement scale, and the CKM matrix elements, including their allowed 
ranges, are the same as in (27j- We take mj/^ — 3.097 GeV, and the J/ip meson distribution amplitudes [HHH, 



x(l — x) 



1 — Aax(l — x) 
x(l — x) 



1 — Aax{\ — x) 



(17) 



with the decay constant /m = 405 MeV [29| . The shape parameter will be varied between a = 0.7 ± 
0.1 to acquire the theoretical uncertainty from these nonperturbative inputs. It turns out that our results 
are not sensitive to the variation of a. The normalization constants N L,t are determined by J dx^> L,t (x) — 
/ J/v ,/(2v/2AQ, giving the central values N L = 9.58 and N l = 10.94. The form factor Ff K (m 2 J/i: ) = 0.62±0.09 
is quoted from the light-cone sum-rule calculation [3(|, to which 15% theoretical uncertainty has been assigned. 

It is observed that Ajj^ K0 and A[jJ^ k0 are of the same order, confirming the importance of the nonfactorizable 
contribution in color-suppressed B meson decays. Moreover, their real parts cancel, such that Ayj^ K0 = 

-^jJ!pK ^^j/t/jK i s a l mos t imaginary. The amplitudes ^4j/'^o lead to the branching ratio, and ASj^x s ano - 
Aj/^Xs from the correction to the decay, 

B(B° -> J/yjK ) = (6.6+^[+ 3 2 ;g) x 10~ 4 , 

As j e /7K S - (7.2^; 4 4 Sf))xio- 4 , 
4mL - -fie.r^S.iOxio- 4 , (18) 



J/^Ks - \- L "-'-8.7(-4.1) 

where the errors in the parentheses arise only from the variation of the hadronic parameters. The result for the 
B° J/ipK° branching ratio is in agreement with the data B(B° -> J/ipK°) = (8.72±0.33) x 10~ 4 [Hj]. A form 
factor F BK (m 2 j^) increased by 15% and larger Gegenbauer coefficients in the kaon distribution amplitudes can 

easily account for the central value of the data. Our ASj e ^ Ks (^j/^Ks) from the penguin operators is about 
twice of the naive estimation from the w-quark-loop contribution [f| , and has an opposite (the same) sign. Both 
AS j e /7K s and A< j/7k s in Eq- TO indicate the O(10" 3 ) penguin pollution in the B — > J/ i^Ks decay, consistent 
with the conjecture made in (l7| . 

We have checked that the correction to the mixing does not modify the denominator 1 + \Xj^x s \ 2 m ^ ne 
definitions of Sj/^k s arj d Aj/^x s m Eq. (jTJ) . Hence, it is legitimate to simply add the values in Eqs. ([2|) and 
Q to Eq. (|18p . We then derive ASj/^k s an d Aj/^x s m the most complete analysis up to leading-power in 
1 /rrib and to next-to- leading order in a s , 

as JNiKs = Asf;; Ks + asj% Ks = (9.3±S:g) x 10- 4 , 

A JHKs = Af;; Ks + AJ% Kg = ~ (14.11?^) x 10- 4 . (19) 

Including the CP violation from the K-K mixing [13], ASj/^^s ano - Aj/^k s remain O(10~ 3 ). Equation (JTUJ) 
will provide an important standard-model reference for verifying new physics from the B° — > J fi^Kg data in 
the future. The estimate of ASj/^k s ~ O(10 -4 ) in [j| implies that they have overlooked the more crucial 
contribution from the penguin operators. Our observation of the direct CP asymmetry supports the claim in 
[3] that Acp(B — > J/tpK) observed at 1% or higher would indicate new physics definitely. 
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